Flavour structure of supersymmetric 5*0(10) GUTs with extended matter sector 
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We discuss in detail the flavour structure of the supersymmetric SO(10) grand unified models 
with the three traditional 16-dimensional matter spinors mixed with a set of extra 10-dimensional 
vector multiplets which can provide the desired sensitivity of the SM matter spectrum to the GUT 
symmetry breakdown at the renormalizable level. We put the qualitative argument that a successful 
fit of the quark and lepton data requires an active participation of more than a single vector matter 
multiplet on a firm, quantitative ground. We find that the strict no-go obtained for the fits of the 
charged sector observables in case of a single active matter 10 is relaxed if a second vector multiplet 
is added to the matter sector and excellent, though non-trivial, fits can be devised. Exploiting 
the unique calculable part of the neutrino mass matrix governed by the SU (2) l triplet in the 54- 
dimensional Higgs multiplet, a pair of genuine predictions of the current setting is identified: a 
non-zero value of the leptonic 1-3 mixing close to the current 90% C.L. limit and a small leptonic 
Dirac CP phase are strongly preferred by all solutions with the global-fit x 2 -values below 50. 

PACS numbers: 12.10.-g, 12.15.Ff, 12.60.Jv, 14.60.Pq 



I. INTRODUCTION 

Even after 35years since the pioneering work by Georgi 
and Glashow [l[ the idea of grand unification still re- 
ceives a lot of attention across the high energy physics 
community, providing one of the most popular schemes 
beyond the Standard Model (SM) of particle interactions. 
Apart from the canonical prediction of the proton insta- 
bility and monopoles, the simplest grand unified theories 
(GUTs) can be tested for the compatibility between the 
observed SM flavour structure and the simplified shape 
of their Yukawa sector emerging at the grand unification 
scale Mq, typically in the ballpark of 10 16 GeV. 

Recently, with the advent of the precision neutrino 
physics [2(, the field experienced a further renaissance 
fuelled by the observation of neutrino flavour oscillations 
3 . The eV scale of the light neutrino masses governing 
these phenomena is often connected to the scale of the 
new physics underpinning a variant of the seesaw mecha- 
nism [4[ . To this end, GUTs can provide a very detailed 
information on the relevant high energy dynamics, with 
implications for the position of the seesaw thresholds and, 
hence, the absolute neutrino mass scale. 

With the new piece of information at hand, the flavour 
structure of the simplest GUTs has been scrutinized thor- 
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oughly in the past [3-0] ■ The intriguing pattern of flavour 
mixing in the lepton sector, together with the constraints 
on the absolute neutrino mass scale, turned out to be ex- 
tremely useful in discriminating among the simplest po- 
tentially realistic GUTs, in particular those based on the 
SO(10) gauge symmetry Q. 

The main virtue of the SO(10) framework consists in 
the fact that every SM matter generation fits perfectly 
into a single 16-dimensional chiral spinor of S'O(IO), thus 
providing a simple rationale for the very special anomaly- 
free pattern of the SM hypercharges. On top of that, the 
right-handed neutrino is inevitable and, hence, seesaw is 
naturally accommodated. As a rank=5 gauge symmetry, 
SO (10) also admits a large number of viable symmetry 
breaking chains resulting in many different interme- 
diate scale scenarios with rich phenomenology. 

From the neutrino perspective, the most important as- 
pect of this freedom is the scale of the B — L symmetry 
breakdown. In the most popular schemes it is triggered 
either by the vacuum expectation values (VEVs) in the 
16-dimensional 15*0(10) spinors or in the irreducible com- 
ponents of the five- index antisymmetric tensor (126©126) 
in the Higgs sector. In supersymmetric (SUSY) scenarios 
with 126# ©126//, the i?-parity of the minimal supersym- 
metric standard model (MSSM) emerges naturally as a 
remnant of the 5O(10) gauge symmetry [13, HH, there 
are no proton-dangerous d = 4 operators and a poten- 
tially realistic Yukawa sector with a calculable seesaw 
can be implemented at the renormalizable level fl2l [l3| . 

On the other hand, one has to resort to a cumber- 
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some Higgs sector as further multiplcts are needed to 
break through an intermediate SU(5) symmetry which is 
left intact by the SM singlets in 126 h © 126//. Remark- 
ably, none of the simplest options, i.e., neither extra 45 
nor 54, is sufficient to do so at the renormalizable level 1 
, and even with both of them, non- renormalizable op- 
erators are still needed to mix the 5[/(2) /.-doublets in 
126// with those from other Higgs multiplets (10- or 120- 
dimcnsional) in order to get a reasonable Yukawa sector. 
Actually, renormalizability calls for 210// instead which 
can provide both the SU(5) breakdown as well as the 
doublet mixing, yet retaining a high level of predictiv- 
ity. Unfortunately, the minimal renormalizable SUSY 
50(10) model [Hf with 10// © 126// © 126// © 210// in 
the Higgs sector does not seem to work due to the generic 
tension between the neutrino mass scale and SUSY uni- 
fication constraints [l5|. Recently, there have been sev- 
eral attempts to overcome this issue by, e.g., invok- 
ing split SUSY [l|| or by employing a 120-dimensional 
Higgs representation (see for instance [l7|, EH and refer- 
ences therein). However, most of these constructions are 
plagued by the instability of the perturbative description 
due to a Landau pole emerging close to the GUT scale 

mm- 

The situation in models with 16// © 16// triggering 
the B — L breakdown [2(| is quite different in sev- 
eral aspects. First, a concise Higgs sector of the form 
16// © 16// © 45// © 54// is sufficient to break through 
the SU(5) lock 2 . Second, here there is no problem with 
mixing the 5£/(2)£-doublets in 10// (which is again in- 
troduced for the sake of a potentially realistic Yukawa 
sector) with those in 16// © 16// at the renormalizable 
level. Moreover, the Landau pole is safely postponed be- 
yond the Planck scale. 

In spite of these attractive features, it turns out to 
be rather difficult to construct predictive models along 
these lines in practise. The basic reason is that there is no 
way to communicate the information about the SU (2)/?© 
U(1)b-l and SU(5) symmetry breaking (driven by the 
VEVs of 16// ©16// and 45 H © 54//) to the matter sector 
spinorial bilinears 16 m 16m at the renormalizable level. 
Thus, in order to get potentially realistic effective quark 
and lepton spectra and mixings, non-renormalizable op- 
erators must be invoked and there is a need for further 
assumptions to retain predictivity in the Yukawa sector, 
see, e.g., @ and references therein. 

An elegant solution [2l| to this conundrum consists 
in abandoning the "matter in spinors" paradigm of the 
50(10) model building. With extra 10-dimensional 
50(10) matter vectors in the game (to be denoted by 



1 In this respect, the situation in the non-supcrsymmetric setting 
differs substantially from the supersymmetric case, c.f. [T3l . 

2 This statement, however, is not a trivial analogue of a similar 
mechanism at play in the 126 © 126 case because the product 
16h54h16/j, unlike 126h54// 126//, docs not contain a gauge 
singlet and thus one of the parameters is missing here. 



10m) admixing at a certain level into the light mat- 
ter fields, the basic invariants of the form l% M lO M l& Hl 
10m10a/54# and 10a/10a'/45// do the magic at the renor- 
malizable level. Moreover, since the SM-singlet VEV of 
16//, (16//), governing the mixing between the spinors 
and vectors can be comparable to the scale of the (gauge 
singlet) mass term MioIOmIOm, the matter vectors do 
not need to decouple from the electroweak-scale (v) 
physics - it's not the v over Miq but the (16//) over Miq 
ratio that matters. 

This is even more so in the SUSY GUTs where a single- 
step breaking (bringing (16//) to the vicinity of the GUT 
scale Mg) is typically favoured. Furthermore, if one ad- 
mits a hierarchy in the eigenvalues of even a Planck- 
scale Mio that could originate from a similar source like, 
e.g., the hierarchy of the Yukawa couplings, it's very 
plausible to expect at least one of them at around (or 
even below) Mq. This, indeed, makes observable non- 
decoupling effects of the extra 10m 's very natural. Re- 
markably, in such 0j C0jSG, the relative magnitude of the 
SU (2)/f © U(1)b-l and 5C/(5) breaking observed in the 
MSSM matter spectra (of the order of the differences in 
the second to third generation mass ratios, i.e., few per- 
cent) is nicely linked to the hierarchy of the SUSY GUT- 
scale thresholds. Moreover, the triplet contribution to 
the neutrino mass matrix turns out to be calculable in 
this framework because the leptonic SU(2)l doublets in 
10a/ can couple to the Higgs triplet in 54// at the renor- 
malizable level. 

Let us also note that the extra vectors in the mat- 
ter sector are inevitable in the unified models beyond 
5O(10), like, e.g., in E 6 GUTs 0. Recently, the ex- 
tra matter in the 5O(10) GUT context played a central 
role in works [23[ in which a class of phenomenologically 
viable models of tree-level gauge mediation as means of 
SUSY breaking has been constructed. 

Although this framework has been used before by sev- 
eral authors to address, e.g., the flavour problem of the 
SM or to constrain the SUSY flavour and CP structure of 
its GUT- inspired extensions (24|, a generic study of the 
flavour structure of the SUSY 5O(10) GUTs with vector 
multiplets in the matter sector has been carried out only 
partially, namely for a single vector matter multiplet at 
play in [25J where a no-go for the simplest setting has 
been formulated. In this study we attempt to go beyond 
the minimal case and look at the viability of a more re- 
alistic scenario in which a hierarchy in the Planck-scale 
Mio brings a pair of its eigenvalues to the vicinity of 
the GUT scale. As we shall see, the generic no-go of 
[25| is lifted already for the second lightest eigenvalue of 
A/10 contributing with just around 1% of the strength of 
the first one and, even within such a "quasi-decoupled" 
setting, the flavour structure of the SM charged matter 
sector is accommodated in a very natural manner. 

Remarkably, complete fits including the triplet- 
dominated neutrino sector observables require a signi- 
ficant contribution from the second 10m hi the matter 
sector, far from the quasi-decoupled regime. In such a 
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case, the minimality of the Higgs potential is fully ex- 
ploited and two generic predictions of the scheme can be 
identified: the best fits of all the measured quark and 
lepton flavour parameters strongly favour small but non- 
zero value of the leptonic reactor mixing angle within 
the ballpark of the current global upper limit [261 ] . to- 
gether with a closc-to-zero value of the leptonic Dirac 
CP phase. 

The work is organized as follows: In section II we define 
the basic framework, derive the effective mass matrices 
for the MSSM matter fields and comment on the role the 
calculable triplet contribution plays in the neutrino mass 
matrix. After a brief recapitulation of the no-go for the 
minimal setting, these formulas are subject to a thorough 
numerical analysis in section III for the case of a pair of 
non-decoupled 10a/ 's and we comment on the blindness 
of the best % 2 fits to the contributions associated to the 
Yukawa coupling of 45// observed in a large part of the 
parametric space available to good charged sector fits. 
In section IV, we briefly comment on the prospects of a 
realistic model building and its basic strategics. Then we 
conclude. 



II. THE FRAMEWORK 

Let us begin with a definition of the minimal frame- 
work in which the generic principles advocated above can 
be implemented in a potentially viable manner. Since the 
details of the matter sector flavour structure depend only 
loosely on the specific shape of the Higgs sector, we shall 
focus on the simplest conceivable model. The following 
discussions can be then extended to more complicated 
settings in a straightforward way. In order to keep the 
discussion compact, we shall stick to salient points only 
and, whenever appropriate, refer to work [25[ where a 
similar construction has been discussed in great detail. 



A. The model definition 

1. The matter sector 

We shall consider the standard three copies of the 
50(10) spinors 16 M (i = 1,2,3) in the matter sec- 
tor (otherwise one could not accommodate properly the 
three generations of up-type quarks), together with n 
copies of the 50(10) vectors 10 M (k = l,..,n). The sub- 
script M indicates that these multiplets are odd under a 
Z^ 1 matter parity invoked in order to prevent the clas- 
sical trouble with the d = 4 proton decay due to their 
potential mixing with the Z^-even Higgs multiplets car- 
rying a generic subscript H. The effective matter sector 
spanned non-trivially over both 16m's and IOm's then 
exhibits a full sensitivity to the GUT-scale VEVs, over- 
coming the "high-energy blindness" of the purely spino- 
rial matter in the renormalizable settings with 16//©16//. 



The SU(3) C © 5/7(2) L ® U(1) Y structure of these mul- 
tiplets reads (in the Q = T| + Y convention): 

16 M = (3, 2, +i)8 (l,2,-i)®(3,l,-|) (1) 

© (3,1,+i)® (1,1, +1)®(1, 1,0) 
10 M = (3,l,-±)ffl(l,2,+±)ffi(3,l,+±)ffl(l,2,-±) 

The SM sub-multiplets of 16m above will be, from now 
on, consecutively called Ql, Ll, E/£, D c L: 7V£ and £"£, 
while those of 10a/ as Al, A l , A c l and A^. 

Let us reiterate that at the SU(3) C <g> SU{2) L ® U{l) Y 
level D L can mix with A C L and Ll with A/, giving rise to 
the physical down quark and chargcd-lcpton components 
(to be called d L and II), sharing the features of both 16a/ 
and 10a/, in particular their sensitivity to the GUT-scale 
physics). 

Let us also note that the matter sector spanned on 
16a/'s and IOm's can be viewed as a hint of an under- 
lying Eq gauge structure where these multiplets both fit 
into its fundamental 27-dimensional representation (de- 
composing under 5O(10) as 27 = 16 8 10 8 1). On the 
other hand, this correspondence is rather loose here as we 
do not demand the number of IOm's to match the num- 
ber of 16m 's, let alone the absence of the extra singlets, 
c.f. section Hi C 11 

2. The Higgs sector 

Concerning the Higgs model that can support the de- 
sired 5O(10) -> SU{3) C © SU(2) L (g) U(1) Y symmetry 
breaking chain at the renormalizable level, the simplest 
such setting in the SUSY context corresponds to the 
16//ffil6//ffi45//©54// Higgs sector. Note that 45// alone 
is not enough because the F-flatncss aligns its VEVs with 
the SM singlets in 16// 8 16// leaving SU(5) unbroken 

m- 

The relevant factors consist of the following SM com- 
ponents: 

16// = (3, 2, +1)8 (1,2, -|) 8(3,l,-f) 

© (3, 1, +|)8(1, 1, +1)8 (1, 1,0) , 

16// - (3, 2, -|) 8 (1,2, +|)8(3,1,+|) 

© (3, 1,-1)® (1,1, -1)9 (1,1,0) , (2) 

45// = (1,3,0)© (1,1, +1)8 (1,1,0) © (1,1,-1) 

© (8,1,0)© (1,1,0) © (3,1, +|)© (3,1, -f) 

© (3, 2, -f ) © (3, 2, +i) © (3, 2, +§) © (3, 2, -±) , 

54// = (1,1,0) ©(1,3,0) ©(1,3, +!)©(!, 3,-1) 

© (6, 1, +|) © (6, 1, -f ) © (8, 1,0) © (3, 2, +±) 

© (3, 2, -|) 8 (3, 2, -i) 8 (3, 2, +|), 

where the underlined SM singlets are all expected to 
receive GUT-scale VEVs. These we shall call V 16 , 
V 16 , V^ 5 (the one in (15,1,1)45 with respect to the 
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SU(4) C ® SU(2) L <g> 517(2)* c 50(10)), V^ 5 (the one in 
(1,1,3)45 i n the same notation) and V 54 , respectively 3 . 
The ultimate SU(2) L ® f7(l)y J7(1)q breakdown is 
then driven by the 5i7(2)z,-doublets in 16// © 16// to- 
gether with a pair of extra copies coming from an addi- 
tional 5O(10)-vector Higgs multiplct 
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H 



(3,1, 



■|)©(l,2,+|)©(3,l,+i)® (1,2,-i), 



which is added as usual in order to end up with a po- 
tentially realistic Yukawa sector. In a self-explanatory 



notation, we shall use the symbols v 



16 

d ■ 



,16 



and u 



in 



for the corresponding doublet VEVs. Apart from these, 
the interplay between the B — L and the electroweak 
breakdown gives rise to a pair of induced VEVs on the 
electrically neutral components of (1,3, ±1) of 54// (to 
be called w±). Subsequently, the renormalizable cou- 
pling 10m10m54# gives rise to a set of Majorana entries 
in the relevant neutrino mass matrix, c.f. section Til Bl 



3. The renormalizable Yukawa superpotential 

The Yukawa superpotential of the model under consid- 
eration reads (with all indices and the Lorentz structure 
suppressed) : 



W Y = 16 m YW h IQm + 16 M 7T6//10 M 
+ 10 M (A54 JJ +t)45h+M 10 )10 m 



(3) 



where Y is a 3 x 3 complex symmetric Yukawa matrix, F 
is its 3 x n general complex analogue in the mixed 16m~ 
10m sector and Mio and A (and rf) are n x n complex 
symmetric (antisymmetric) matrices. At the 5t/(3) c (g> 
U(1)q level, the part of our interest can be written as: 



,in 



W Y 3 U L YUtv^ + NlYN L vH 
+ D L YDlvf+E c L YE L v™ 
+ D L FA c L vf + E c L FA~ L vf + NtFAfv™ 



(4) 



+ D C L FA L V 16 + E L FA C +V ie + N L F AfV 16 
+ A° L \A° L w + +Af\Afw- 
- A L XA C L V 54 + f A^+AA-V 54 + f Af AA° V 54 
+ A LV A C L V 45 + A C L + V A- L V£> + A? V A° L Vf 
+ A L M 10 A C L + A c +M w Al + A?M ia A° L , 

where the defining SU(2)l doublets have been broken 
into their components, i.e., Ql = {Ul,Dl), Ll = 



3 Recall that V 16 is connected to V 16 by the desired D-flatness of 
the SUSY vacuum: |V 16 | = |V 16 |. The notation for the singlets 
in 45// is justified by the observation that V^ 5 can give masses 
only to the quark-like states in 10jv/' s while V^ 5 enters only the 
leptonic bilincars. This is clear from the Pati-Salam decompo- 
sition of the 5O(10) vector which reads 10 = (6, 1, 1) © (1, 2, 2) 
where the former factor a 
corresponds to A/ © A5 . 



(N L ,E L ), A L = (A° L ,Al) and A C L = (A c + ,Af). Wher- 
ever possible, we have also absorbed the relevant Clebsch- 
Gordan coefficients into 0(1) redefinitions of the inde- 
pendent VEVs and/or couplings, with an important ex- 
ception at line 6 where the ratio of the Clebsches can not 
be hidden. This, indeed, is the backdoor through which 
the desired SU (5) symmetry breaking due to a non-zero 
V 54 is transferred into the matter sector. 



B. GUT-scale mass matrices 

The relevant GUT-scale mass matrices for the matter 
fields can be readily read out of eq. (0]): 





= Yv* , 




(Yvf 


M d 


\f t v 16 


M e 


(Yv? 





Fvf 



/0 Yv™ 


v- • 




(5) 
(0) 

(7) 
(8) 



For the first three (Dirac) mass matrices above, the 
following bases have been used: {Ul){Ul) for M u , 
(D L: A L )(Dl,AD for M d and (E L , A£ )(££, A c +) for 
M e , respectively. The Majorana mass matrix M„ has 
been given in the symmetric basis (Nl, iV£, A°, A^°). We 
have also made use of the symmetry properties of Y, Mio, 
A and rj and defined 



Ma 



Mio - AF 54 + r,V£ 5 , 
Mio + |AV 54 - 77 V A 45 



(9) 



Inspecting the matrices above one can appreciate the role 
of the extra vector multiplets in propagating the informa- 
tion about the intermediate symmetry breaking into the 
matter sector: First, since there are no heavy partners 
to the up-type quarks available the physical spectrum 
is determined solely by the spinorial bilinear Yukawa Y. 
Second, the hierarchy of the down-type quark spectrum is 
clearly different from the up-type quarks whenever there 
is a non-negligible admixture of the A C L components in 
the light eigenstates. For this to be the case, | J^V 16 1 
should not be negligible with respect to Ma- Third, in 
order to account for the differences in the down-quark 
and charged-lepton mass hierarchies it is inevitable to 
have Ma different from Mj which can happen only if at 
least one of the S'C/(5)-breaking VEVs V£>, V£ 5 and/or 
V 54 is turned on and it is not screened by the 5O(10)- 
singlct mass term Mio m ©• Thus, at least some eigen- 
values of Mio are required to be in the vicinity of the 
GUT scale. Note that, in spite of the 5O(10)-singlet 
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nature of Mio, this can easily be the case if Mio hap- 
pens to exhibit a several-orders-of-magnitude hierarchy 
as some other Yukawa couplings in the game, in particu- 
lar Y oc M u . 

Note also that there are several interesting formal lim- 
its in which the matter spectrum reveals an enhanced 
symmetry pattern: 

• Putting V^ 5 = V^ 5 and V" 54 to zero one has M A = 
Ml and thus M d = AlJ due to the residual SU(5) 
gauge symmetry left unbroken by V 16 . 

• For Mio strongly dominating the heavy sector 
masses, the extra vectors 10^ decouple and the 
sensitivity of the light sector to the intermediate 
symmetry breaking scales is lost. In this case, all 
Dirac masses are proportional to each other due to 
the residual SU(4)c Pati-Salam symmetry exhib- 
ited by the matter sector, as expected in all settings 
with 16mY16m10h alone in the Yukawa sector. 

• For V 16 < M a ,a with M 10 , U 54 and U 45 at around 
the GUT scale the effect of the SU(A)c symme- 
try breaking becomes observable only in the heavy 
sector because of the effective suppression of the 
FV 16 term linking the GUT-scalc VEVs to the 
light cigenstates. In other words, the vector matter 
does again decouple from the 50(10) spinors. 

These remarks demonstrate clearly the internal consis- 
tency of formulas (O-©. 



C. Effective mass matrices 

Below the GUT scale the heavy part of the matter 
spectrum decouples and one is left with the three stan- 
dard MSSM families. Their masses and mixings are then 
dictated by their projections onto the defining basis com- 
ponents 16 M and 10 M , providing the desired sensitivity 
to the GUT symmetry breakdown in the matter sector. 

In what follows, we shall use the calligraphic symbols 
M f (with f = u,d, e, v) for the effective MSSM mass 
matrices to make a clear distinction between these and 
the full- featured GUT- level mass matrices ([S])-©. 



1. Integrating out the heavy degrees of freedom 



a. Up-type quarks: Since there are no multiplets in 
the 10m, c.f. decompositions (TTJ), with the up-type quark 
quantum numbers the effective MSSM up-quark mass 
matrix (evaluated at the GUT scale) is identical to the 
5O(10)-level mass matrix ([5]): 



M u 



Yv, 



10 



(10) 



b. Down-type quarks and charged leptons: The situ- 
ation is very different though for down-type quarks and 
charged leptons whose GUT-level mass matrices ^ and 
(O are (3+n) x (3+n)-dimensional. They can be brought 
into a convenient form by means of transformations 



M d 



M' d , M e 



U* e M e = M' e , (11) 



where U d , e are (3 + n) x (3 + n) unitary matrices such 
that M' d and M' e are block-triangular: 



ML = O 



v v 
Mq 



, M' e = Q 



v 

V Mr 



(12) 



This corresponds to the change of basis in the right- 
handed (RH) down quark and left-handed (LH) lepton 
sectors respectively: 




(13) 



Here the upper components of the rotated vectors (d c L 
and II) correspond to the light MSSM degrees of free- 
dom. Note also that the residual SU (2) l gauge symme- 
try makes the GUT-scale rotations (fT3")) act on both the 
charged lepton (El; A£) as well as the neutrino (Nl\ 
A L ) components of the leptonic doublets Ll and A^. 

Since the residual rotations acting on the LH quark and 
RH charged lepton components bringing the M' d e matri- 
ces into fully block-diagonal forms are extremely tiny (of 
the v/Mq order of magnitude) the 3x3 upper-left blocks 
(ULB) in relations (fT2j) can be readily identified with the 
effective light down-type quark and charged lepton mass 
matrices, i.e., M d = (A^)ulb, M e = (MQtjlb- Given 
the specific form of M d and M e in eqs. (O and ([7]) and 
parametrizing the unitary matrices Ud and U e as 



U d . e = 



C d . K D 



d.c 



d,e ^d,e 



(14) 



(here A d<e , B de , Cd. e and D de are 3x3, 3xn, nx3 and 
n x n matrices, respectively) one obtains 



YAlvf+FB\vf, 



M T e = YA\vf + FBlvf. 



(15) 
(16) 



The off-diagonal GUT-scale blocks of M d and M e are 
rotated away provided 



F T A d V 16 + M A B\ = 0, 
F T AlV w +MlB\ = 0, 



(17) 
(18) 



which link the A dfi and B d e factors. The last two rela- 
tions, together with the unitarity of Ud t e implying 



A dte A\ e + B d , e Bl e = l, 



(19) 



impose strong constrains on the elements of matrices (TTJ} 
entering the effective mass formulas (jT5"T) and ([TBI . These 
correlations shall be fully exploited in section IIIII 
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c. Neutrinos: The situation in the neutrino sector is 
slightly more complicated due to the higher dimension- 
ality of the GUT-level mass matrix ([SJ. Notice, however, 
that the action of the LH leptonic rotation ([13")) , corre- 
sponding to a transformation M v — > U*M V U], = M' u 
with U v denoting the relevant (6 + 2n) <g> (6 + 2n)- 
dimensional unitary matrix, yields Ml in a hierarchical 
form 4 



M' = 



(B* e XBtw + A* e Yvl° B* e \Dtw + 
YC\v™ 

d;xdIw + 

V • 



Fvf 
M K 
Xw-J 



, (20) 



with an abbreviation Mj = F T C\V W + M\D\ for the 
only GUT-scale entry therein. 

Naively, given the hierarchies of the SU(2)l triplet, 
doublet and singlet VEVs, this shape of M' v yields three 
electroweak-scale pseudo-Dirac neutrinos at the effective 
theory level (corresponding to the upper- left 6x6 block of 
M' v above), in an obvious conflict with observation. This 
is namely due to the fact that the lower-right (3 + 2n) x 
(3 + 2n)-dimensional sector of matrix (|20| corresponding 
to the (7V£,A°,A^°) part of the rotated basis does not 
have a full rank at the GUT scale. 



2. Calculable triplet seesaw 

However, this issue should not be taken very seriously 
unless the quantum stability of the small entries in M' v is 
discussed. In particular, the 22 block zero (correspond- 
ing to the -ZV£7V£ bilinear in Wy, i.e., a SM singlet-singlet 
contraction) is not protected by the electroweak symme- 
try and thus can be naturally subject to large corrections 
which, eventually, may restore the full (GUT-scale) rank 
of the lower-right block of M' v . 

For instance, a dimension 5 operator of the form 
16a/ 16Ml6.HT6.fr/Mp, where Mp is the Planck scale, lifts 
this zero sufficiently to change the entire picture: the 
lower-right block becomes superheavy and the hierarchi- 
cal matrix structure a la standard seesaw is achieved. 
Subsequently, one is left with three sub-eV Majorana 
neutrinos at the SM level, with the upper-left entry of 
M' v promoted to the role of an additive (type-II-like) con- 
tribution to their effective mass matrix. 

Let us also remark that a simple renormalizable re- 
alization of this scheme is obtained if the matter sec- 
tor is further extended by three 5*0(10) singlets, well in 
the spirit of Eq gauge models. The extra contraction 
16m1m16jj in the Yukawa superpotential provides the 
necessary set of large matrix elements entering the heavy 



part of the (extended) neutrino mass matrix even at the 
renormalizable level. 

However, given the likely proximity of such a new 
physics scale to M G , one expects other physical ef- 
fects to affect all the effective mass matrices at some 
level. Obviously, it is not very appealing to let the non- 
renormalizable operators and/or similar effects into play 
in the simple scheme of our interest unless these are under 
a very good control 5 . Actually, as we have already em- 
phasized, the goals of the current analysis are rather dif- 
ferent and, as long as we focus on the renormalizable part 
of the effective flavour structure, a deep understanding 
of all the neutrino sector details is not strictly required. 

Indeed, whatever the ultimate rank-restoration mech- 
anism happens to be, the seesaw contribution due to the 
SU(2) L triplet in 54 H , 



M^ = B* e XB\w + 



(21) 



is always present and the underlying 10m10a/54# con- 
traction is particularly robust. Indeed, apart from the 
standard SU{2)l <£> U(l)y gauge symmetry protection, 
this is namely due to the fact that the triplet VEV within 
54ff can not be mimicked by (45^) nor (16^-) at the d = 5 
level. Since A4„ is also the only calculable part of the 
effective neutrino mass matrix in the simple framework 
of our interest, the best one can do is to focus entirely on 
it and assume its dominance over the other contributions 
in M v : 



Mv~M*. 



(22) 



This approximation is what we shall adopt from now on. 
Let us also note that a dedicated analysis of the con- 
ditions under which such situation can be realized in a 
specific complete model is a highly non-trivial enterprise, 
much beyond the scope of this work. 



III. ANALYSIS AND DISCUSSION 

A. General prerequisites and comments 

With all this information at hand one can attempt to 
exploit the strong correlations between the effective mass 
matrices (flU)) , (|T5|). (fTo]) and (|2"T) to assess the viability 
of the general framework by means of a global % 2 analysis 
of its compatibility with the measured quark and lepton 



masses and mixings. 



4 Note that the upper-right corner zero is due to the SU(2)l gauge 
symmetry which promotes the requirement l|18|l of a similar zero 
in the charged lepton mass matrix mil l to neutrinos. 



Note that giving up renormalizability one would actually loose 
a great deal of the original motivation for the vector-like matter 
entering the genesis of the SM flavour structure, as discussed in 
section [I] Indeed, there is a lot of non-renormalizable SO(10) 
models of flavour in the literature with spinorial matter only. 
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1. The effective quark and lepton mass matrices 

Before that, one should attempt to further simplify the 
relevant mass matrices <jT2J) - (HSJ) • First, one can substi- 
tute Y for Ai u and eliminate the Bd, e factors in eqs. 
(|T5|) . (|T6)l and (21]) by using relations dTTJ) and CEH]) so 
that AdgS remain the only "complicated" factors in all 
formulas of our interest: 

M d = (rM u — F(A~I A )~ 1 F T ^J A a , (23) 

Ml = (rM u - F(Ml)- l F T ) A\, (24) 

M u oc A* e F(M A )- 1 ~X(Ml)- 1 F T Al, (25) 

with r = vf/v™, M u ,d,e = M u ,d, e /m b , F = 
Fvf/m h , Ma, a = {vf /m b )M AA /V w and A = 
{vf/m b ){V 54 /V 16 )\ where m b stands for the bottom 
quark mass. In what follows, it will also be convenient to 
normalize the antisymmetric parts of Ma, a in the same 
manner: t^a.a = (v^ 6 /mb)(V A 5 A /V 16 )ri. Note that the 

overall scale of M. v driven by w+ remains undetermined 
at the current level. For this reason we have dropped 
the explicit triplet VEV and introduced a proportional- 
ity sign into eq. (|2"5"|) . 

It should be also possible to write down the Ad, e fac- 
tors in terms of the Yukawa supcrpotcntial parameters 
as we did for the brackets in eqs. ([23]) and (|24[) which, 
however, could be quite complicated in general . Actu- 
ally, we don't need to do so as there are redundancies 
in Ad >e that do not play any role in the low energy phe- 
nomenology (i.e., spectra and LH mixings). Indeed, one 
can always decompose Ad, e as 

A d , e = Vd^Hd.e (26) 

where Vd. e and Hd. e arc unique 3x3 unitary and hcr- 
mitcan matrices respectively. It is clear that Vd, e do not 
affect the low energy quark and lepton observables be- 
cause Vd contributes only to the RH quark rotations and 
V e enters Ai e and M. v on the same footing and thus can- 
cels in the lcptonic mixing matrix. 

Given (|2!?|) the hcrmitean factors Hd,e can be deter- 
mined from the unitarity of Ud and U e (|19j) taking into 
account the triangularization constraints (IT7|) and (fT5|) : 6 

H d = (1 + F*(M A M A )~ 1 F T )~ 1 / 2 , (27) 
H e = (1 + F*{MlM* A )- 1 F T )- 1/2 . (28) 

To conclude, formulas (|2"3l) - (j2"51) admit for a full recon- 
struction of the quark and lepton masses and mixing (up 



The square root of a generic hermitian positive scmidefinite ma- 
trix M is defined as U\/~Dlj' { where D = U^MU is a real non- 
negative diagonal matrix. Note that the sign ambiguity in \[D 
does not play any role due to the irrelevance of the overall signs 
of the generalized eigenvalues of matrices I)23|l - lj25|l and the cor- 
responding mixing angles. 



to the absolute neutrino mass scale and irrelevant ba- 
sis transformations V d , e ) for any point in the parametric 
space of the model. 

2. Basic features and strategy for potentially realistic fits 

Let us now comment on the salient features of the ef- 
fective flavour structure ([23 ]) -([25 ]) and its prospects for 
accommodating successfully the quark and lepton data. 

• First, it is clear that for non-zero F and Ma. a the 
up and down quark mass matrices as well as the 
hierarchies of their spectra arc different and a non- 
trivial quark mixing is generated. 

• The Cabibbo-Kobayashi-Maskawa (CKM) quark 
mixing angles are naturally generated when the 
magnitude of the F(M A )- 1 F T term in (J33]) is 
smaller than that of rA4 u ] otherwise the approx- 
imate alignment of M. u and Ai d is lost and there 
is no reason for the CKM mixing to be small. 

• In such settings, the r parameter has a clear in- 
terpretation of a "hierarchy compensator" between 
rrib and rrit and as such its value is strongly con- 
strained. As we shall recapitulate in section UlI B 21 
(c.f. [25j), pushing r out of this natural domain 
hampers the prospects of getting good fits of both 
the down quarks and the charged leptons at once. 

• The case of sub-leading F{M A )~ 1 F T and 
F(Mj)" 1 F t naturally accommodates the ap- 
proximate convergence of the b and r Yukawa 
couplings observed in many studies of the running 
of Yukawa couplings. 

• Moreover, for F in the 0(1) ballpark, the same 
implies F*(M A M A )- 1 F T ,F*{MlMl)- 1 F T < 1, 
which provides a further insight into the effective 
mass formulas for Aid and Ai e because it renders 
the A d and A\ factors in (|23|) - ([24|) unimportant 
even for the second generation. 

• With ^(Ma)" 1 ^ and ^(Mj)- 1 ^ 71 in a few per- 
cent domain there should be enough room to ac- 
commodate the differences among m c /m t , m s /mb 
and m^/m T . Moreover, even the basic hierarchy 
between the CKM mixing angles #12 ^ #23. 13 seems 
very natural: with a diagonal Ai u the only CKM 
angle that can be large due to a few-percent off- 
diagonalities from the sub- leading term is #12- 

• The neutrino mass matrix (|25[) has nothing to do 
with the leading contribution to Ai e and thus there 
is no reason for the leptonic mixings to be small. 

Remarkably, this scheme matches perfectly the basic 
qualitative features of the observed quark and lepton 
mass and mixing pattern. In what follows, we shall be 
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using the values given in TABLE Q] as physical inputs of 
the numerical analysis carried out in section fill CI 



Quark sector 


observable 


value 


observable 


value 


m u [MeV] 


0.45(±0.2) 


m d [MeV] 


1.3 ±0.6 


m c [MeV] 


217(±35) 


m a [MeV] 


23 ±6 


mt [GeV] 


97(±38) 


m b [GeV] 


1.4 ±0.6 


sin 6\ 2 


0.2243 ± 0.0016 


sin 0*3 


0.0351 ± 0.0013 


S in0? 3 


0.0032 ± 0.0005 


°CP 


60° ± 14° 


Lepton sector 


observable 


value 


observable 


value 


Amlj [eV 2 ] 


(7.7 ±0.2) 10" 5 


m e [MeV] 


0.3565 ± 0.0100 


\Am 2 31 \ [eV 2 ] 


(2.40±0.12)10 -3 


m M [MeV] 


75.3 ± 1.2 


sin 2 e[ 2 


0.304 ±0.019 


m-r [GeV] 


1.629 ±0.037 


sin 2 6 l 23 


0.50 ±0.06 


sin6»i 3 


< 0.18 



TABLE I: Sample GUT-scale inputs of the numerical analysis 
performed in sectio nTH CI The specific values correspond to 
those given in [53] for the quark sector and [28|] (c.f. also 
[29l ]) for the charged lepton masses, tan/3 = 55°. The solar 
and atmospheric neutrino mass squared differences and the 
leptonic mixings are taken from The upper bound on 6[ 3 
corresponds to the global 90% C.L. value quoted in [2(|. The 
running effects in the neutrino sector have been neglected due 
to the hierarchical shape of the neutrino spectrum. For sake 
of simplicity, symmetric cr-ranges have been adopted. The 
error in the electron mass has been artificially enhanced by a 
factor of 10 to improve the convergence of the numerics, with 
no significant impact on the quality of the actual fits. 



3. Parameter counting 

In order to assess the prospects of testing this picture 
even at the quantitative level it is worth counting the 
number of independent parameters. Working with a real 
and positive M u (that fixes entirely the basis in the space 
of 5*0(10) matter spinors) the phase of r can be rotated 
away by a global phase redefinition of Aid and A4 e , leav- 
ing a single real parameter (RP). A similar rotation in 
the space of n 5O(10) matter vectors can bring the M w 
matrix to the real and diagonal form with n RPs. In this 
basis, the complex symmetric Yukawa coupling of 54/r 
(A) adds (n ± l)n RPs and the antisymmetric Yukawa of 
45# (i]), which is present for n>2, yields (n — l)n RPs. 
For n > 2, one must also add the complex ratio of the 
two VEVs in 45//, accounting for an extra pair of RPs. 
Finally, there is the 3 x n-dimcnsional complex matrix of 
F's adding in general 6n RPs. In total, one ends up with 
2n 2 ± 7n + 3 RPs for n > 2 (and 10 RPs for n = 1, in 
agreement with (25j). 

With the up-quark masses as inputs, there are 13 low 
energy observables one can attempt to fit (3 down-quark 
masses plus 4 CKM parameters in the quark sector, 3 
charged lepton masses, the Am^/IAm^l ratio in the 



neutrino sector and 2 leptonic mixing angles measured so 
far). A successful fit of these data could then admit to tell 
something about the unknown parameters (in particular, 
sin B\ 3 and the leptonic CP phases). 

B. Single active vector matter multiplet 

Let us begin with the case of a single vector matter 
multiplet in the game 7 . For n = 1, however, the neutrino 
mass matrix (|25[) has rank 1 and thus there is no point 
in attempting to fit Amfi/IAml-J nor 6 l l2 . Hence, in full 
generality, one is left with 10 parameters to fit 11 observ- 
ables, which clearly indicates a potential difficulty with 
a full-fledged three-generation fit. Nevertheless, since in 
practice there can easily be other 10m 's around (though 
perhaps at the verge of decoupling) it still makes sense 
to look at the two heavy generations. As we shall see in 
section fill B II an interesting link between the maximal- 
ly of the atmospheric mixing in the lepton sector and 
the interplay among the 23 mixing in the quark sector 
and the m s /mb ratio can emerge even in this obviously 
oversimplified case. Moreover, in order to appreciate the 
naturalness of the n = 2 fits discussed in section IIII CI 
it is instructive to see explicitly where the trouble with 
the three-generation fit for n = 1 [25[ comes from; an 
analytic argument will be given in section IIII B 21 

1. Triplet seesaw and a large 2-3 mixing in the 2x2 case 

Perhaps the most intriguing feature of the minimal sce- 
nario is the simple correlation between the large values of 
the leptonic 2-3 mixing inherent to the triplet-dominated 
neutrino masses and the specific flavour structure ob- 
served in the 2-3 part of the quark sector. It reads 

tan 20^ w 2\x\j |l- x 2 \ , (29) 

with x = (y b /y s ) sm $23' where y St b are the Yukawa 
couplings of the heavy down-type quarks (in the di- 
agonal basis) and sin 0| 3 is the 2-3 mixing angle in 
the quark sector. In a certain sense, this relation 
can be viewed as a "radiatively stable" analogue of 
the well-known Bajc-Senjanovic-Vissani (BSV) relation 
tan 26^ 3 w sin 20f 3 /2 sin 2 Q\ z ± e (with e = 1 - y T /y b ) [H 
derived in the minimal SUSY 5O(10) GUT framework 8 . 
The relation of our interest (f2T?)) is readily obtained from 
the basic formulas for the charged lepton and the triplet 
neutrino masses (|2"4"|) and (|25[) taking into account the 
estimated structure of the charged sector fits specified in 



7 Since this case has been analysed in detail in [TjEll here we shall 
just recapitulate the salient features of this basic setting. 

8 What we mean here by "radiative stability" is that the yb/ys 
ratio is subject to a much milder running than the ratio y T /yb 
underpinning the BSV relation in the minimal SUSY SO(10). 
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section UlI A 21 or in j25|. For sake of simplicity, we shall 
also assume a CP-conserving setting with all phases ei- 
ther or 7r. At the leading order, the flavour structure of 
the triplet-dominated neutrino mass matrix can be ap- 
proximated by 

M u oc B* e \B\ oc A* e FF T A\ « V*FF T V^ , (30) 

where wc made use of the fact that A is a number now 
and the "external" factors A e are almost unitary, see sec- 
tion [nrX2] Rotating away the V e matrices, the charged 
lcpton mass matrix (|24|) becomes close to diagonal. Thus, 
focusing entirely on the 2-3 mixing (which, indeed, is 
the only leptonic angle it makes sense to look at with a 
rank=l mass matrix), it is almost entirely encoded in the 
neutrino mass matrix M. v oc FF T and one can write 



relations for A 



for these factors: 



tan 2023 ~ 2l.F2.F3l 



Fi 



(31) 



In order to get a grip on the typical values of the F- 
paramcters in (|31[) one should exploit the quark sector 
sum-rule. At the same level of accuracy as before, the 
relevant formula (f2"3"|). once contracted to the 2nd and 
3rd generations, yields 



M d 



m c jvn,b 






m t /m b 



F 2 F2F3 
F2 F3 Fo 



(32) 



where p = (m& / 'v^ 6 )V 16 j 'M a and the approximate diag- 
onality of Aid in the .A/f^-diagonal basis has been used. 
Due to the estimated smallness of r (in the few % range), 
one can expect (c.f. section IIII A2[) that the only rele- 
vant entry of the first matrix in (|32j) is rat /nib. The 
resulting Aid can be easily shown to give m s /mb ~ PF2, 
1 m rmt/rrib and sin 023 ~ PF2F3. Solving for F2 and F3 
and substituting into (|3Tt one recovers (|29|) . 



2. The renormalizable 3x3 charged sector no-go 

With such an observation at hand one would naturally 
ask whether the analysis can be extended to the 3x3 
case so that it might account for the details associated to 
the light flavours. Unfortunately, the answer is negative. 
The reason is that with a single vector matter multiplct 
at play there is a fundamental obstacle to any poten- 
tially successful fit already at the charged sector level. 
Remarkably enough, one can even provide a simple ana- 
lytic argument for why this happens to be so. 

For the sake of that, let us look at the shape of the 
down-quark mass matrix (|23j) and consider the three 
main minors of AtdA4 d defined as Aj 



da djj 



di-jdji = dudjj - \dij\ 2 with dij = (M d M d )ij- Notice 
that these quantities, by definition, depend only on the 
physical inputs, in particular the quark masses and mix- 
ing parameters. One can easily show that the F couplings 
enter A^^j only as |Fi| 2 |F.,-| 2 (recall there is only a sin- 
gle 10m here so F is a vector) and one can solve the three 



p 2 \F t \ 2 



FA Z = 



{da 



r 2 mF] 



j j 



r 2 m? 2 ) 



+ m?u - 27774 mi cos(7 l - 7j) 



(33) 

where ji are the phases of Fi defined as Fi = \Fi\e 'w 2 . 
It is clear that consistency requires the numerators on 
the RHS of cq. (f3"3")l to be non- negative for all i, j. First, 
this can never be realised nontrivially if the CKM mix- 
ing was turned off (implying dy^^y = 0) - at least one 
pair out of any three non-zero numbers always yields a 
positive product. Thus, with a single 10m at hand, a 
non-trivial Vckm is a necessary condition for any suc- 
cessful quark sector fit. Second, turning on the small 
CKM mixing, the numerators look like products of pairs 
of quadratic functions in r 2 with small positive shifts due 
to \dij\ 2 0. Taking into account the physical ranges of 
the quark masses and mixing angles it is straightforward 
to check that the only domain, in which all three of these 
expressions can be simultaneously positive, corresponds 
to r k m s /m c ~ 0.15. This value, however, is one order 
of magnitude away from the physically motivated expec- 
tation identified in section IIII A 21 at odds with the de- 
sired shape of the charged lepton spectrum. Thus, even 
at the pure charged sector level, there is a generic no-go 
for the fits of the flavour structure of the minimal model 
with a single vector matter multiplet [25j . 

On the practical side, one should emphasize that the 
argument above is based on the specific values of the 
input parameters used throughout this analysis, c.f. TA- 
BLE HI These, however, depend on a particular scenario 
employed to study their running properties. For instance, 
large SUSY thresholds [3l[ can significantly affect the 
GUT-scale mass ratios, especially for the light genera- 
tions, and, hence, the desired range for the r parameter. 
Thus, at least in principle, there could still be an option 
for the n = 1 case to be implemented in models yielding 
unconventional high-scale Yukawa patterns. A detailed 
discussion of these issues, however, is beyond the scope 
of this work. 



C. Two active vector matter multiplets 

In view of the negative result for 77 = 1 it is natural to 
ask whether the charged-sector no-go can be overcome 
with more than a single extra matter multiplet in the 
game, in particular with 77 = 2, and if yes how much one 
can learn about the leptonic mixing (especially about 9 l 13 ) 
and CP violation in such case. At first glance, one would 
expect the 77 = 2 fits to be essentially trivial as the dimen- 
sionality of the parametric space increases dramatically: 
from 10 for 77 = 1 to 25 for 77 = 2, c.f. section llll A 31 On 
the other hand, the 3 extra constraints from the leptons 
can play an important role, given the qualitative differ- 
ence among the hierarchies and mixings in the quark and 
lepton sectors. Moreover, as we know from the previous 
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section, good fits are impossible if the second 10m plays 
only a marginal role, i.e., if it dynamically decouples. 

To put this statement on a firm ground one should 
take into account how the gauge-singlet mass parameters 
encoded in the M 10 matrix enter the heavy matter spec- 
trum. In the normalization A/io = (v^ 6 /m^AIio/V 16 
(see sectio nTlI Al) one can conveniently parametrize 



M w = tdiag(l,p) 



(34) 



where t is an overall factor and p encodes the hierarchy of 
the two eigenvalues of M^. Note that one can take p > 1 
without loss of generality. Then, p — > oo corresponds to 
the decoupling limit if the couplings between the heavy 
and the light GUT-scalc matter states arc kept undcr 
control, i.e., do not diverge. The expected worsening of 
the best \ 2 towards the decoupling limit can then be 
used as a non-trivial consistency check of the numerical 
results we shall present in the subsequent sections. 



1. Fits with 45 h decoupled from the Yukawa sector 

Let us begin with the case of a negligible contribution 
from the 45/j Yukawa coupling. As we shall see in sec- 
tion MI C 2| this is well motivated because of a high de- 
gree of "sterility" of 45// in the n = 2 fits whenever there 
is more than an 0(1) hierarchy between the mass terms 
of the two lOjvf- Note also that the situation with rj — s- 
is effectively parametrized by only 21 RPs; thus, taking 
into account the strong phenomenology constraints on r 
together with the perturbativity bounds on l-Fifcl's and 
the limited impact of their phases, it is actually far from 
clear whether this setting admits good fits. In what fol- 
lows, we shall use a simple prescription for the relevant 
decoupling parameter 9 , 



P = p/max{|c 22 |, |ci 2 | 2 } , 



(35) 



where cm govern the entries of the properly normalized 
Yukawa coupling of 54# , 



Afej = t Cki 



(36) 



The shape of formula (|3"5j) reflects the basic features of 
the numerical fits, namely the dominance of |c22 1 <x P for 
large p followed by a milder behaviour of |ci2 1 oc ^fp and 
an essentially p- insensitive \c\\ | cx p°. Apart from the ci 2 
playing the obvious "destructive" role of mixing up the 
heavy and the light sectors, the c 22 is taken into account 
because it can mimic an "effective" p in Ma or M\. 



9 The specific form of the definition Il35jl corresponds to the role the 
C22 and C12 factors play in the spectrum of a which (for cu 
in the 0(1) domain) is well approximated by their determinants 
(linear in p and C22 and quadratic in C12). 



-J 



iog 10 (P) 



FIG. 1: An n = 2 histogram of the relative frequency of 
fits obtained for the quark and charged lepton masses and 
CKM mixing parameter with < 1 for different values of 
the decoupling parameter P. The sharp decline of the counts 
towards the high P limit is a manifestation of the no-go dis- 
cussed in section IHIB 21 for the n = 1 case. The shaded 
region on the left corresponds to the fine-tuned setting with 
V 54 dominating over the singlet mass parameter Mio and its 
specific shape is an artefact of the numerical method we use. 



a. Pure charged sector fits - avoiding the n = 1 no-go: 
The first test to be passed concerns the charged lepton 

fits in the n = 2 case. Recall that in section UlI B 21 these 
were shown to be generally troublesome in the n = 1 case 
despite the relatively large number of parameters (10 in 
general) available to fit just 7 observables (3 down-type 
quark masses and 4 CKM mixing parameters if, for the 
sake of simplicity, the up-type quark masses are fixed at 
their means). With the extra 10m at hand, excellent fits 
are easily obtained within the expected domains (see sec- 
tion MI A 2p whenever its contribution is non-negligible. 
Quantitatively, as seen in FIG. [TJ we have found good 
fits of the charged sector data for all values of the decou- 
plings parameter P below about one hundred. In other 
words, the value of P^ 1 ~ 1% constitutes a qualitative 
boundary above which the second vector matter multi- 
plet is already decoupled too much to avoid the no-go 
inherent to the n = 1 settings. 

In this respect, it is also worth noting that the inter- 
esting link between the large atmospheric mixing and the 
specific value of mb jm s sin #23 sa 1 obtained in the n — 1 
case, c.f. section MI B 11 is upset due to the perturba- 
tions coming from the second vector matter multiplct 
and there is no preferred value of 8 l 23 observed in these 
fits. 

b. Fits including leptonic 9\ 2 , # 23 and Am^i/IAm^l: 
Including from now on the measured values of the rel- 
evant neutrino oscillation parameters, i.e., 8 l 12 , #23 ano - 
Am^/IAm^l, into the \ 2 function one can still attempt 
to get predictions for 8[ 3 and the leptonic Dirac (S l CP ) 
and Majorana CP phases. Remarkably enough, such fits 
turn out to be nontrivial in spite of the high number 
of free parameters at play. This is reflected by the fact 
that none of the fits we obtained yields \ 2 below around 
15; nevertheless, given the number of fitted observables, 
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these values are still to be regarded as very good. 

The generic behaviour of the relevant fits can be seen 
in FIG. [2] Now, the best \ 2 value is a steeply rising 
function of the decoupling parameter P, in agreement 
with expectation. On the more technical side, here we 
have also decided to lift some of the residual degenera- 
cies in the parametric space by further constraining the 
Cu parameters into the 0(1) domain, which provides a 
convenient link between the p and P parameters. A de- 
tailed information about a pair of the relevant best-x 2 
solutions is given in TABLE UU 




FIG. 2: A sample of the \ 2 values for the n — 2 fits with 
77 — > as a function of the decoupling parameter P. One can 
see clearly that the extra constraints from the leptonic sector 
make the n — 2 fits non-trivial, see also section MI C 21 and 
TABLE HU 



2. Complete fits including 45h 

Turning on the antisymmetric Yukawa coupling of the 
one could expect that the extra parameters associ- 
ated to this sector would make the global fits of the mea- 
sured quark and lepton masses and mixing parameters 
much simpler than in the 77 — > case discussed above. 
On the other hand, it is also clear that 45// should have 
almost no impact in the quasi-decoupled regime when, 
effectively, only one of the two IOm's contributes to the 
light states. Thus, the situation is more subtle and, as 
one can see in FIG.[3j the extra 45//-term in the mass for- 
mulas actually leads to a significant improvement of the 
fits only in the very-low- P region of the parametric space, 
c.f. FIG. [3J As before, a detailed information about a 
pair of low-% 2 solutions can be found in TABLE IIIII In- 
deed, since P and p are again strongly correlated, c.f. 
section MI C II in both cases also the p parameter falls 
into the 0(1) domain. 

a. Sterility of 45// for p > 10: 
Although the tight link between p and P emerging in the 
|cjy| ~ 0(1) regime justifies the high degree of sterility 
of the 45// contribution for large p ~ P > 0(100) values 
corresponding to a quasi-decoupling of the second 10m, 
it could be rather surprising that very good fits can be 



Parameter 


Fit I 


Fit II 


P 


14.339010 


2.847552 


r 


0.0f621f50 


0.01473598 


t 


31.791794 


162.846941 


rh u 


0.0003270355 


0.0003340457 


rh c 


0.f6f8762 


0.1698742 


rh t 


69.008012 


78.470173 


Fu 


-0. f 69686 - 0.f63037i 


-0.398016 - 0.499151? 


F21 


0.568262 + 1.543220i 


-0.322017 - 0.584528i 


F31 


-2.992396 + 1.508342i 


4.433113 - 1.105430? 


F12 


0.300966 + 0.871563i 


0.997104 - 0. 219207i 


F22 


-1.198678 - 2.f97259i 


3.814132 + 0.370594i 


F32 


-1.731386 - 0.f58647i 


7.384473 - 3.331683i 


Cll 


1.958581 - 3.831928i 


-0.766690 - 1.509670i 


C22 


-0.550307 - 0.871218?! 


-0.031916 + 0.033495i 


C12 


2.555502 - 3.497512i 


0.769329 - 1.704622? 


m u [MeV] 


0.4579 


0.4677 


m c [MeV] 


226.6 


237.8 


mt [GeV] 


96.61 


109.86 


m d [MeV] 


0.8892 


0.9909 


m a [MeV] 


40.24 


30.50 


m b [GeV] 


1.461 


1.634 


sin 9\ 2 


0.2248 


0.2240 


sin 0f 3 


0.03487 


0.03153 


sin 9\ 3 


0.003304 


0.003958 




37.38° 


60.83° 


m e [MeV] 


0.3561 


0.3582 


m M [MeV] 


75.29 


75.32 


m T [GeV] 


1.630 


1.588 


|Am^ | 


0.03269 


0.03244 


sin 2 d[ 2 


0.2714 


0.3031 


sin 2 e l 23 


0.3323 


0.4207 


2 

Xtotal 


21.319 


15.222 



TABLE II: A sample pair of low-^ 2 solutions in the n — 2 
case with 77 — ► (c.f. section fill C f bp . The four digit accu- 
racy adopted in the physical parameters reflects the maximum 
quality of the input data these quantities are compared to, c.f. 
TABLE HI Let us also remark that a full reconstruction of the 
displayed \ 2 values an interested reader could attempt could 
be partly obscured by the limited precision of the displayed 
numbers. 



obtained only for p ~ 0(1). A thorough inspection of 
the role of 45jy in the relevant mass formulas given in 
Appendix [X] reveals that this is namely due to the anti- 
symmetry of the corresponding Yukawa coupling 77 which 
gives rise to, e.g., a further 0(m s /uib) suppression of the 
45// effects in some of the quark sector obscrvables, in 
particular the first and second generation masses and the 
f 3 and 23 CKM mixing angles. 
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FIG. 3: A sample of the x 2 values for the full n — 2 fits as a 
function of the decoupling parameter P. One can see clearly 
that the extra constraints from the leptonic sector make the 
n = 2 fits troublesome, see also section [III C 21 



b. Genuine predictions for 9 l 13 and 8 l CP : 
In the fits above, we let only the well measured quark and 
lepton masses and mixing parameters contribute to the 
global x 2 function. The other observables, in particular 
the reactor mixing angle and the leptonic CP phases were 
left apart as genuine predictions of the current scheme. 
Indeed, for any specific fit, these can be calculated in 
terms of the other parameters listed in TABLE IIIII 

In FIG. [4] we display the predicted values of the lep- 
tonic 13 mixing obtained for the fits indicated in FIG. [3] 
with x 2 iS 150. Although it is impossible for the best-x 2 
points to get within the current 90% C.L. experimental 
limit, there is a clear preference of a small 1-3 mixing 
with sin Si 3 ~ 0.2 at the low-x 2 tail of the distribution. 



1 p 

0.8- 
0.6- 
0.4 
0.2=.* 

0- 



C 

03 



..... ... ..." ' i. JT* IJug&iStr* 



50 



100 



150 



FIG. 4: The predicted value of the leptonic 13 mixing as a 
function of the x 2 corresponding to the fits of all the other 
measured parameters. The current 90% C.L. upper limit 
sin #5.3 < 0.18 is indicated by the dashed line. The distri- 
bution of the calculated 9 l 13 values for the lowest-;^ 2 points 
clusters in the lower part of the available domain at around 



Parameter 


Fit I 


Fit II 


P 


3.041675 


2.847552 


r 


0.01312141 


0.01499823 


t 


87.744176 


162.609350 


rh u 


0.0003216498 


0.0003551049 


rh c 


0.1461787 


0.1823098 


rht 


80.619489 


77.955444 


#ii 


0.006594 - 0.01261H 


-0.397570 - 0.498592i 


#21 


1.868096 + 0.406222i 


-0.320084 - 0.581020i 


#31 


-0.865185 + 12.270027i 


4.510621 - 1.124757i 


#12 


0.192324 + 0.789105i 


0.976804 - 0.214744? 


#22 


2.488878 - 0.413259i 


3.817823 + 0.370953i 


#32 


0.321369 + 2.679685i 


7.491695 - 3.380059i 


Cll 


-7.019216 + 4.756008i 


-0.766062 - 1.508434; 


C22 


-2.349768 - 0.927223i 


-0.032584 + 0.034197; 


C12 


0.2841363 - 0.567570z 


0.769877 - 1.705836; 


dA 


-0.988269 + 0.783068i 


0.121933 + 0.169316i 


d A 


0.005013 + 0.001892i 


0.015713 + 0.000770; 


m u [MeV] 


0.4503 


0.4971 


m c [MeV] 


204.7 


255.2 


m t [GeV] 


112.9 


109.1 


m d [MeV] 


0.6364 


1.0176 


m s [MeV] 


26.97 


30.79 


m b [GeV] 


1.268 


1.651 


sin 9f 2 


0.2241 


0.2248 


sin 6»| 3 


0.03465 


0.03289 


sin 6\ z 


0.003243 


0.003546 


°CP 


47.77° 


55.13° 


rn e [MeV] 


0.3562 


0.3571 


7Ti M [MeV] 


75.30 


75.30 


m T [OeVJ 


1 cm 

t.6ty 






0.03226 


0.03226 


sin 2 e[ 2 


0.3039 


0.3048 


sm 2 e l 23 


0.4769 


0.4152 


2 

Xtotal 


3.203 


9.187 


sin 9[ 3 


0.269 


0.255 


S'cp 


-10.57° 


11.99° 



sin 



0.2. 



TABLE III: A sample pair of low-x 2 solutions obtained 
in section IIII C 21 for the n = 2 case. Here d/\,A = 
f _1 (v^ 6 /mi,)(V r ^ 5 A /V rl6 )?7i2 and the last two rows represent 
the relevant predictions for 8[ 3 and 5 l CP obtained with the 
corresponding fits, c.f. also FIG. [4] and FIG. [5] 



Similarly, as one can see in FIG. [SJ a small leptonic 
Dirac CP phase S l CP is strongly preferred for the lowest- 
X 2 solutions. As far as the Majorana phase is concerned 
(recall that one of the light neutrinos is exactly masslcss 
in the current setting) we do not observe any specific fea- 
ture in its distribution and the predictions are essentially 
uniformly covering the whole available domain. 



D. More than two vector matter multiplets? 

As we have seen, even with 2 copies of extra matter 
multiplets at play the fits of the system (f2"3"|) - ([25]) are 
non-trivial, although a naive parameter counting (see sec- 
tion IIII A 3[) would clearly suggests the opposite. As a 
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FIG. 5: The predicted value of the leptonic Dirac CP phase 
5 l CP (in units of 7r) as a function of the \ 2 corresponding to the 
fits of all the other measured parameters. The distribution of 
the calculated 5 l CP values for the lowest-;^ 2 points clusters at 
around zero, thus indicating a possible difficulty in revealing 
the leptonic CP violation in the next generation of neutrino 
oscillation experiments. 



matter of fact, this is very welcome because such a set- 
ting admits to draw genuine prediction that can be tested 
at near-future experimental facilities. 

From an underlying Eq perspective one could ask 
whether a third 10m would cause a qualitative change 
of the picture. Given the number of extra parameters 
popping up in the n = 3 setting the general answer is 
very likely to be positive. For the same reason, this is 
not the strategy we would like to pursue as it would most 
probably lack any predictive power. Moreover, reiterat- 
ing the hierarchy arguments given in section HI the over- 
all scale of the corresponding Mio would be unnaturally 
low if one brought all three of its eigenvalues below the 
Planck scale. Apart from these rather technical issues, a 
third 10m would not shed any new light onto the neu- 
trino sector which, in spite of its appeal, doesn't need to 
be dominated by the triplet contribution at all. 

Hence, without a handle on the non-renormalizablc 
terms governing the type-I sector, we consider further 
(n > 3) extensions of the current analysis to be rather 
academic. Nevertheless, there are various concepts that 
can provide an extra information making such studies 
nontrivial and potentially interesting, be it family sym- 
metries, extra dimensions, finite unifications or anything 
else. This, however, is beyond the scope of this study. 



[23| where it was shown to be capable of accommodat- 
ing a simple tree-level realisation of the gauge-mediated 
supersymmetry breaking mechanism in a phenomenolog- 
ically viable manner. 

Focusing on the next-to-minimal case with more than 
a single such matter 10 playing an active role in the effec- 
tive SM matter spectrum, the well-known no-go emerging 
in the minimal model already at the charged sector level 
is alleviated. Subsequently, we have pursued an exten- 
sive x 2 analysis including the neutrino sector data which 
is tractable only if the calculable 5t/(2)L-triplet contri- 
bution dominates the neutrino mass matrix. We have 
obtained very good fits of all quark and lepton masses 
and mixing parameters measured so far, providing a pair 
of genuine predictions for those to be, presumably, within 
the reach of the near future facilities: the reactor mix- 
ing angle is predicted to be relatively larg e, close to the 
current 90% C.L. limit quoted, e.g., in (2q |. while the lep- 
tonic Dirac CP phase tends to be very small and, hence, 
more difficult to access. No preference has been observed 
for the relevant Majorana CP phase. 

Unfortunately, a full account of the neutrino sector in- 
cluding also the type-I-like contribution to the seesaw for- 
mula is intractable without an additional information on 
how the full rank of the heavy part of the neutrino mass 
matrix is restored in a specific setting. For instance, one 
can think about extra flavour symmetries 32[ that may, 
at least to some extent, keep the number of free param- 
eters under control, and at the same time constrain the 
Yukawa couplings of the model. This, however, typically 
requires non-renormalizable operators to be invoked at 
some level, thus challenging the original motivation for 
the extra vector-like matter as a renormalizable key to 
the observed quark and lepton masses and mixing. Nev- 
ertheless, at closer look such models are likely to exhibit 
rather different correlation patterns due to the generic 
dominance of better-controlled renormalizable contribu- 
tions. In this respect, an extra flavour symmetry could be 
at least partially unloaded from the usual burden of hav- 
ing to address many aspects of the SM flavour problem 
at once, as it is often required when matter is spanned 
over the 5*0(10) spinors only. 



IV. CONCLUSIONS AND OUTLOOK 
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Appendix A: Sterility of 45h for p > 10 



As seen in section llll C 1 b| in the 77 — >• limit the good 
fits of all observables including the neutrino ones required 
a very mild hierarchy between the two eigenvalues of M\q, 
i.e., p < 10. Remarkably, such a strong preference of very 
low p values appears also in the general fits with non-zero 
77 where one would expect it to be much weaker due to 
the extra freedom associated to the active role of (45//) 
in the relevant formulas. 

Although it is quite difficult to provide a general un- 
derstanding of this behaviour one can address at least 
some of its aspects. In particular, one can decipher why 
for a given moderate-x 2 point obtained in the 54j/-only 
fits of section IIII C II the extra freedom associated to the 
subsequent inclusion of 45/r in section IIII C 21 does not 
improve the \ 2 of the corresponding complete fits. 

First, it is easy to show that for 2x2 matrices the anti- 
symmetric part of the inverse of an arbitrary matrix M 
is a function of the antisymmetric part of M only, apart 
from an overall normalization. Note that this specific 
to 2 x 2 matrices and does not hold for larger dimen- 
sionality. Thus, the antisymmetric fj enters the formulas 
(f23]) - (|24p in a very specific manner: it only generates an 
extra antisymmetric contribution to the ubiquitous sym- 
metric part of the FM~ 1 F T bilinears generated by the 
A/10- an d A- pieces in (properly normalized) eqs. ©. As 
such, barring the sub- leading effect it has in Hd and H e , 
it can affect the relevant mixing angles at the 0(e) level 
while the spectrum of the matrices (|23 |) - ([24|) remains in- 
tact up to 0(e 2 ) corrections [TtJ where e is parametrizing 
the "smallness" of the antisymmetric correction as com- 
pared to the symmetric one given by the remaining terms 
in 1(23 ]) - ([24"| ) . Since the FM^ A F T bilinears are tailored 
to give rise to the second generation masses, the typi- 
cal size of their leading order entries in the down-quark 
sector is m s /mi, while it is m^/m T for the charged lep- 
tons, both in the few percent range. However, this is all 
namely due to its symmetric part dominated by M\q and 
for p > 1 there is an extra overall suppression associated 
to the antisymmetric piece. This can be seen, e.g., from 



(Ma)- 1 = 



detM A 




(Al) 



where s k i = (M10 - \)kl, a = 1)12 and detM A = sns 2 2 - 
s 2 2 + a 2 . It is clear that for moderate A and fj of the order 
of 0(t) the leading contribution to the symmetric part 
of FM~ 1 F T scales as \F\ 2 s kk /s 11 s 2 2 « \F\ 2 /t while the 
antisymmetric piece \F\ 2 a/sns 2 2 ~ \F\ 2 /pt is suppressed 
by an extra factor of p^ 1 . Thus, the relevant e parameter 
behaves like m s jm^p ~ 0.02/p. From this, it is already 
clear that for the fits with p > 10 such an antisymmetric 
correction can not help lowering the \ 2 value of a specific 
fit if it comes predominantly from the second generation 
of down quark or charged lepton masses. 

Concerning the impact of a non-negligible 77- 
contribution to the mixing parameters the situation is 
somewhat more subtle. At the leading order, one can 
quantify the shift in the CKM mixing angles as: 



Vckm 



V, 



CKM 



(1-Z). 



(A2) 



where Z is an anti-hcrmitcan matrix obtained from 



Aij/{Sd)j 



and A stands for the antisymmetric part of F(Ma) 1 F t 
in the basis in which the symmetric part of Md, Sd ~ 

{rM.u — F{M^)~ 1 F T y is diagonal and real. Since, as 
we have seen, A is l/p-suppressed with respect to the 
symmetric part of F(Ma) _1 -F t , one has 



Z ps - 



/0 0(1) 0(<5) N 
. 0(6) 
\. . 



with 5 — m s /m& , (A3) 



where we have approximated the eigenvalues of Sd by 
{md/rrib, m s /mb, 1}. Thus, for the "transition region" 
values of the p-parameter, i.e., p > 10 under consider- 
ation, only 6\ 2 can be slightly affected by the Yukawa 
of 45/f while the other CKM mixings remain essentially 
intact. This, however, does not improve the fits with 
"moderate" x 2 value which is spanning over several dif- 
ferent observables other than just 0\ 2 . 

Let us also remark that the situation in the lep- 
tonic sector is very similar to quarks, in particular for 
the charged lepton contribution to the leptonic mixing. 
Moreover, the current precision of the leptonic mixing 
parameters determination is much worse than the same 
in the quark sector so the net effect of the antisymmetric 
Yukawa of 45# in the leptonic mixing x 2 contribution is 
essentially negligible for p > 10. 
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